Getting the astrophysics and particle physics of dark matter out of next-generation 

direct detection experiments 
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The next decade will bring massive new data sets from experiments of the direct detection of 
weakly interacting massive particle (WIMP) dark matter. Mapping the data sets to the particle- 
physics properties of dark matter is complicated not only by the considerable uncertainties in the 
dark-matter model, but by its poorly constrained local distribution function (the "astrophysics" of 
dark matter). I propose a shift in how to think about direct-detection data analysis. I show that by 
treating the astrophysical and particle-physics uncertainties of dark matter on equal footing, and by 
incorporating a combination of data sets into the analysis, one may recover both the particle physics 
and astrophysics of dark matter. Not only does such an approach yield more accurate estimates of 
dark-matter properties, but it may illuminate how dark matter coevolves with galaxies. 
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Dark matter makes up ~ 23% of the energy density of 
the observable Universe, yet its identity is unknown (e.g., 
[H). While microlensing constraints limit the fraction of 
dark matter that consists of compact objects such as pri- 
mordial black holes [2j , there is a plethora of viable dark- 
matter particle candidates, most arising naturally in ex- 
tensions to the standard model (SM) of particle physics 
1- 

If the dark-matter particle is massive (m x > 1 GeV) 
and interacts with baryons, it may be observed in 
direct-detection experiments *4|. This class of particle 
is generically called a weakly-interacting massive par- 
ticle (WIMP), and the experiments look for the small 
(~ 10 — 100 keV) kinetic-energy transfer from a WIMP 
to a nucleus. Current limits on the WIMP-nucleon spin- 
independent and WIMP-neutron spin-dependent cross 
sections come from ~ 10 kg-sized experiments 5, 6]. Ton- 
scale experiments are expected by ~ 2015, and there are 
plans for ~ 10-ton experiments by 2020 [7|. These ex- 
periments are touted as being able to probe far into su- 
persymmetric and Kaluza-Klein parameter space. 

The event rate in these experiments not only depends 
on the particle-physics properties (m x , cross sections) of 
the dark matter, but on the local dark-matter distribu- 
tion function (DF). The event rate (per kg of detector 
mass in which the target nucleus has atomic number N) 
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where d<jN /dQ is the differential scattering cross section, 
/(x,v) is the local WIMP DF, 
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is the minimum speed required for a WIMP of mass m x 
to deposit energy Q to a nucleus of mass m^r, and /j,n 
is the particle-nucleus reduced mass. Thus, any limits 
or constraints on the WIMP particle-physics properties 
derived from the data depend on the WIMP DF. 



However, the local WIMP DF is poorly understood. 
Even if the local WIMP DF is dominated by a smooth, 
equilibrium dark-matter halo, and the density profile 
could be determined (e.g., Ref. |8j), the velocity distri- 
bution cannot be uniquely determined [§[. Dark- matter 
experiments are the only probe of the local WIMP ve- 
locity distribution. Moreover, the smooth halo may not 
dominate the direct-detection event rates. Recent work 
has suggested that the Milky Way may also have a disk 
of dark matter, and it may dominate the direct-detection 
signal [10] . However, its properties depend sensitively 
on the accretion history of the Milky Way, which is un- 
known. In addition, the WIMP velocity distribution may 
be dominated by a few tidal streams. Simulations can- 
not currently resolve such small-scale structures. With- 
out better understanding of the astrophysical proper- 
ties of the local dark-matter distribution, particle-physics 
constraints on dark matter from direct-detection experi- 
ments are essentially meaningless. 

While there has been some attention to how different 
DFs change exclusion curves in the case of no signal, 
there has been little work on how to determine WIMP 
particle-physics properties if at least one direct-detection 
experiment sees a WIMP signal [1, The bulk of 



that work has the underlying assumption that the direct- 
detection signal is dominated by the smooth halo. 

New ideas and methods are required to properly ex- 
tract particle-physics implications from upcoming data 
sets since it is not feasible to model the astrophysical 
properties of WIMPs. I propose treating the astrophys- 
ical uncertainties in the WIMP DF on par with the 
particle-physics uncertainties, in order to derive robust 
constraints on both the WIMP DF and particle physics 
from the data. However, there are degeneracies among 
the particle-physics and astrophysics properties. One 
must analyze multiple data sets together to break the 
degeneracies. This approach is inspired by cosmological 
data analysis (e.g., [l|). 

As a test of principle of how both particle-physics and 
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astrophysics properties may be derived by a combined 
analysis of multiple data sets if WIMPs are detected, I 
perform a simple toy study employing Fisher-matrix fore- 
casting for toy models of upcoming direct-detection ex- 
periments. Future work will use Monte Carlo simulations 
for forecasting instead of Fisher matrices, but Fisher ma- 
trices suffice for this first look into how well both the 
particle-physics and astrophysics properties of WIMPs 
can be inferred from future data sets. For this study, the 
local WIMP DF is treated as a Maxwellian with a one- 
dimensional velocity dispersion w rms , with the Solar Sys- 
tem lagging the WIMP distribution by a velocity Vi ag . I 
place no restrictions on w rms or i>i ag , in other words, there 
are flat priors on the velocity parameters. Unlike Ref. 
[H, there is no prior that the direct-detection event rate 
must be dominated by smooth halo WIMPs. I only as- 
sume that for whatever WIMP structure dominates the 
event rates, the velocity distribution is described by a 
single Maxwellian distribution. I use this simple form of 
the velocity distribution without priors on the parame- 
ters because this work is a test of principle, to show the 
power of a joint analysis without the theoretical prior on 
which WIMP population dominates the event rates. In 
later work I will input more complex velocity distribu- 
tions and extract the velocity distribution in a nonpara- 
metric way, perhaps akin to Ref. (12| . 

The free parameters are m x ; D = p x o^ j(2m^), the 
normalization of the differential event rate, where p x 
is the local dark-matter density, and a^ 1 is the spin- 
independent WIMP-proton cross section; i>i ag ; and v Tms . 

I set p x — 0.3 GeV cm~ 3 for concreteness, although the 
exact value of the WIMP density is not relevant for this 
analysis, as p x is degenerate with a^ 1 . I set the spin- 
dependent cross sections to zero, a restriction I will loosen 
in future work. I use toy models for a set of direct- 
detection experiments with data sets in 2015 and 2020. 
These experiments can distinguish nuclear-recoil events 
from backgrounds on an event-by-event basis, but they 
cannot determine the recoil direction (although dark- 
matter time projection chambers can 13]). The 2015-era 
set of experiments I use are as follows: (1) A xenon-based 
experiment inspired by XENON1T [3] with a nuclear re- 
coil analysis window 2 — 27 keV and ~ 1/6 efficiency of 
finding nuclear recoil events in the 1-ton experimental 
volume over the course of one year. (2) A xenon-based 
experiment inspired by LUX with a xenon mass of 300 
kg Q. I assume a three-year exposure with an analysis 
window 5 — 27 keV and a 1/6 efficiency. (3) A germa- 
nium experiment modeled on SuperCDMS Phase B. I 
assume a one-year exposure, 145-kg target mass, and an 
efficiency and analysis window similar to that of CDMS- 

II (10 - 100 keV) @, El. (4) An argon experiment in- 
spired by WArP, with a target mass of 140 kg, threshold 
of 40 keV, three- year exposure, and efficiency similar to 
that of the prototype. The larger version of WArP is 
different in form from the prototype, so the efficiency 



and analysis window of the actual WArP experiment is 
likely to be different than what I assume here 16]. The 
2020-era experiments used are as follows: (1) A 20-ton 
xenon experiment, whose characteristics I model like the 
LUX-inspired experiment. (2) A 1-ton germanium ex- 
periment with characteristics of the 2015-era germanium 
experiment. (3) A 10-ton argon experiment. All 2020-era 
experiments are assumed to have one-year exposures. 
I perform the analysis assuming m x = 100 GeV, 
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(D = 1.4 x 10~ 4; 



t'lag = 220 km s , and u rms = 155 km s 
lipses are based on Poisson likelihoods, 
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where N e is the expected number of events given the 
theoretical parameters 9, N a is the observed number of 
events with observed energies {Ql,Q„, ...,Q^°), and 
is the likelihood that a recoil i with observed energy Q l a 
is found given the theoretical parameters. The likelihood 
for each individual event is 
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where £ is the efficiency of identifying an event with 
energy Q l a . For recoils outside the analysis window, 
£(Qo) = 0- For now, I neglect energy errors. The analy- 
sis of the germanium experiments with the CDMS energy 
errors described in Ref. Q is virtually indistinguishable 
from the case of perfect energy resolution. 

In Fig. [TJ I show la errors in the m x -D plane. In the 
left-hand plot, I show errors assuming flat priors in all pa- 
rameters, meaning that there is no external information 
on parameters from other experiments or observations. 
The lightly filled ellipse corresponds to the 2015-era ex- 
periments, and the dark ellipse corresponds to the 2020- 
era experiments. In Fig. [21 the left-hand plot shows errors 
in the v\ ag -v lnls plane with flat priors. While the 2015- 
era experiments do not constrain m x , D, or v\ ag , v rTas 
can be determined to ^50%. The 2020-era experiments 
can constrain all four parameters to of order 20%. 

I compare these errors against those found either with 
a mass prior or assuming that the local WIMP DF is 
perfectly known. In the upper central and right-hand 
panels of Fig. 1, I show error ellipses assuming that the 
LHC has measured m x to 10%, which is the expected 
uncertainty in m x if the LHC identifies a supersymmet- 
ric WIMP 17]. All parameters but m x have flat priors. 
While there are degeneracies among parameters in each 
individual experiment, the experiments complement each 
other to drive down the uncertainties, as does the mass 
prior. In the lower central and right-hand panels in Fig. 
1, I show the errors assuming that the local WIMP DF 
is known exactly. The errors are deceptively small in the 
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FIG. 1: la error ellipses for m x and D if m x = 100 GeV, 
(jp / = 10 -44 cm 2 , and p x = 0.3 GeV cm -3 . Upper left panel: 
Forecasted errors for 2015-era (yellow or lightest grey) and 
2020-era (black) experiments with flat priors on all parame- 
ters. Upper center panel: Error ellipses assuming the LHC 
has measured the WIMP mass to 10%, for 2015-era exper- 
iments. The ellipses represents errors obtained with the 
WArP-inspired (magenta; background medium grey), LUX- 
inspired (blue; dark grey), XENONlT-inspired (red; fore- 
ground medium grey), and SuperCDMS-inspired (green; light 
grey) experiments. The yellow (lightest grey) ellipse shows 
errors for the joint analysis of all experiments with the 
LHC mass prior. Upper right panel: Error ellipses with the 
LHC mass prior for 2020-era experiments (blue or dark grey: 
xenon, green or light grey: germanium experiment, and ma- 
genta or medium grey: argon). The diagonal line shows the 
errors obtained in a joint analysis of all experiments with the 
LHC mass prior. Lower panels: Error ellipses if the WIMP 
DF is fixed. 



2015 forecast, but the errors are not much smaller in the 
2020-era experiments. However, in both fixed 
local WIMP DF will yield significant biases in m x and 
D as well as too-small error bars 0. 

In Fig. [2j I show la uncertainties in the velocity pa- 
rameters. As in Fig. [I] I show errors for each experiment 
analyzed with the LHC mass prior for the 2015- and 2020- 
era experiments, as well as the errors if all experiments 
are analyzed together with the mass prior. The broad 
degeneracy between v\ ag and v Tms in each experiment is 
broken since the direction of the degeneracy is different in 
each experiment. In the joint analysis, the uncertainties 
in the velocity parameters are only slightly reduced with 
the LHC mass prior than if m x were a free parameter. 

In Figs. [3] and 21 I show la uncertainties for dif- 
ferent, "streamlike" velocity parameters but with the 
same particle-physics parameters as in Figs. Q] and [2] 
If 

v iag — v rms — 30 km s , most nuclear recoils are 
low in energy, so experiments such as XENON IT or 
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FIG. 2: la errors in the velocity parameters. Left panel: Er- 
rors obtained in a joint analysis of all 2015-era (yellow or 
lightest grey) and 2020-era (black) experiments with flat pri- 
ors on all parameters. Center panel: Error ellipses for each 
2015-era experiment analyzed with an LHC mass prior. The 
yellow (lightest grey) ellipse shows the errors for the joint 
analysis. Right panel: Errors for the 2020-era experiments 
and an LHC mass prior. Experiments the same as for Fig. [T] 

LUX, which have low energy thresholds and heavy tar- 
get nuclei, are necessary to characterize the velocity and 
particle-physics parameters, although high-threshold ex- 
periments are not constraining. In Fig. G2 the 2015- 
era experiments do a better job of constraining parame- 
ters than the 2020-era experiments because of the lower 
threshold of XENON1T, although constraints on D are 
still weak even with the LHC mass prior. However, the 
mass prior significantly improves constraints in v Tlas . 

In Fig. 21 I show constraints assuming wi ag = 
400 km s _1 and v lms — 30 km s _1 . In this case, the typ- 
ical nuclear recoil is large [e.g., Eq. (j2J)], so experiments 
with large analysis windows (typical of germanium and 
argon experiments) are best able to constrain parame- 
ters. Without mass priors, the degeneracy between m x 
and D is large for 2015-era experiments, but the veloc- 
ity parameters are better constrained. For the 2020-era 
experiments, the uncertainties on the velocity parame- 
ters are tiny even without mass priors. Parameters con- 
straints improve significantly with the LHC mass prior. 

There are three takeaway points from this initial study. 
(1) It is possible to glean both particle-physics and as- 
trophysics inferences about dark matter by analyzing ex- 
periments together if a signal is seen in at least one. (2) 
It is only possible to do this by doing a joint analysis 
of the experiments. The degeneracies in each experiment 
are broken by the other experiments. Thus, having many 
different dark-matter experiments is a necessary condi- 
tion for being able to extract both particle-physics and 
astrophysics properties of dark matter from data. (3) 
Dark-matter experiments are the only probes of the lo- 
cal WIMP DF. As we come to understand in simulations 
what drives the evolution of the WIMP DF, we can learn 
something about the evolutionary history of the Milky 
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FIG. 3: Same particle-physics parameters as Figs. [T] and [2] 
but with i>iag = Dims = 30 km s _1 . Ellipse colors have the 
same meanings as before. 
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FIG. 4: Same particle-physics parameters as Figs. [T] and [2] 
but with vi ag = 400 km s _1 and v TrnB = 30 km s _ . Ellipse 
colors have the same meanings as before. 



Way from the local WIMP DF as unveiled by these ex- 
periments. 

The joint analysis of dark-matter data sets can and 
should be used to extract the most information possi- 
ble about both the particle physics and astrophysics of 
dark matter from direct-detection experiments and other 
types of experiments and observations sensitive to dark 



matter. This is particularly important because the data 
sets expected in the next decade will be vast and diverse. 
By using robust joint analysis methods, we will be able 
to explore the complementarity of ongoing and upcoming 
experimental and observational programs to determining 
the nature of dark matter. 
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